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The description of many dynamical problems like the particle motion in higher dimensional spher- 
ically and axially symmetric space-times is reduced to the inversion of hyperelliptic integrals of all 
three kinds. The result of the inversion is defined locally, using the algebro-geometric techniques of 
the standard Jacobi inversion problem and the foregoing restriction to the ^-divisor. For a represen- 
tation of the hyperelliptic functions the Klein- Weierstrafi multivariable a-function is introduced. It 
is shown that all parameters needed for the calculations like period matrices and abelian images of 
O ' branch points can be expressed in terms of the periods of holomorphic differentials and ^-constants. 

The cases of genus two, three and four are considered in detail. The method is exemplified by the 
■ particle motion associated with genus one elliptic and genus three hyperelliptic curves. Applica- 

tions are for instance solutions to the geodesic equations in the space-times of static, spherically 
, symmetric Hofava-Lifshitz black holes. 
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I. INTRODUCTION 
A. The mathematical problem 

Various problems of physics are reduced to the inversion of a hyperelhptic integral. Namely, let 

y2 ^ 4^2g+l ^ y^^^^2g + . . . + Aq (I.l) 

be a hyperelhptic curve Xg of genus g with one branch point at infinity realized as a two-sheeted covering over the 
extended complex plane. A point P S Xg has coordinates P — (a;, y), where the sign of the second coordinate y 
indicates the chosen sheet on a Riemann surface. Let TZ{x, y) be a rational function of its arguments x and y. We 
consider here the problem of the inversion of an abelian integral 



f TZ{x,y)dx^t (1.2) 



resulting in a function x{t) which is a function of the complex variable t. The integral (jl.2p can be decomposed by 
routine algebraic operations to 



9 pX 9 pX nX 

£{x) - £{xo) + ^ak / duk + ^bk / dr^. + ^ / dfl^^ 

k = l "'^0 k=l "'^0 fe=l "'^0 



t , (1.3) 



where £{x) is an elementary function including logarithms and rational functions, ak, bk, Ck are certain constants and 
dui^jdr^. and dfi^^^^ are differentials of the first, second and third kind, respectively. Namely, duk are holomorphic 
differentials, dr^ are meromorphic differentials of the second kind with a unique pole of the order 2g — 2k + 2, and 
d^a^^p^. are meromorphic differentials of the third kind with first order poles in the points ak and /3k and residues 
±1 in the poles. We suppose that in the case considered there are n > Q differentials of the third kind. 

It is well known that only in the case of elliptic curves, i.e. for g — I, the correspondence x ^ t is one-to-one and 
the aforementioned inversion problem can be solved in terms of single-valued elliptic functions. In the case of higher 
genera, g > 1, a one-to-one correspondence is achieved between the symmetrized products of curves Xg x . . . x Xg 
and a multi-dimensional complex space, the Jacobi variety. Single-valued functions in this case appear to be multi- 
periodic functions of many complex variables, called abelian functions. These ideas already developed by Jacobi had 
led Riemann to the concept of the Riemann surfaces, to the introduction of multi-dimensional ^-functions, and the 
formulation of his celebrated theorems. 

Moreover it is also known that the function x{t) becomes single-valued on the infinitely-sheeted Riemann surface. 
In particular, when abelian integrals reduce to elliptic integrals, this Riemann surface becomes finitely-sheeted [GSOTi] . 
For the special case of a genus two hyperelhptic curve a deta iled construction of su ch an infitely-sheeted Riemann 
surface surface was given by Fedorov and Gomes-Ulate FG07 |. It also follows from |FG07| that x{t) is well defined 



on the complex plane from which an infinite lattice of polygons with (4g — 4)-edges called "windows" is extracted. In 
our investigation we are considering a function x(t) defined on the complement to these windows and suppose that 
the integration paths in (jl.3p never intersect these windows. 

Our approach to the inversion of hyperelhptic integrals is based on the well developed theory of hyperelhptic 
abelian functions and on various relations between the ^-functions and ^-constants. We are describing the function 
x(t) as the restriction of an abelian function, which can be expressed in terms of symmetric functions of the divisor 
in the associated Jacobi inversion problem, to the one-dimensional stratum of the 0-divisor. We are implementing the 
Klein- Weierstrafi realization of hyperelhptic functions in terms of multi-variative cr-functions that represent a natural 
generalization of the standard Weierstrafi g-function to hyperelhptic curves of higher genera. 

This paper continues our recent work |EHKKL"Ti | where the inversion of hyperelhptic holomorphic integrals was 



considered. The novelty of our approach is the simultaneous consideration of all three kinds of abelian integrals from 
a unified viewpoint. At the heart of our method lie algebraic expressions of the symmetric bi-differential of the second 
kind involving the explicitly given Kleinian 2-polar. The inversion procedure involves the theta-constant relations (the 
Thomae and Bolza formulae) that manifest the link between branch points and Riemann period matrices. We are 
also developing a computer algebra procedure that allows to use Maple/algcurves software without explicit knowledge 
of a homology basis encrypted in the Tretkoff-Tretkoff algorithm. This will enable us to find all needed quantities 
like the vector of Riemann constants, the period matrices of the first and second kind, as well as the correspondence 
between branch points and ^-characteristics. We emphasize that in this investigation we concentrate on the algebraic 
side of the derivation The function x(t) obtained in that w ay is defined only locally while its analytic continuation 
never intersects the infinite set of cuts introduced in |FG07( . 
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Another approach to the problem of inversion of integ rals of t he second and third kind that is ba sed on the 
generalized ^-function goes back to Clebsch and Gordan CG866j | and was developed in [Fedod lBF08l| . Here we 
do not discuss generalized Jacobians, and we plan to make a comparison between these two methods in another 
publication. 



B. Physical motivation 



The mathematical results described in this paper have direct applications to the solution of the geodesic equation 
in certain Hofava-Lifshitz black hole space-times. The Hofava-Lifshitz theory [H09al lH09bl | is an alternative gravity 
theory that is powercountable renormalizable. The basic idea is that only higher spatial derivative terms are added, 
while higher temporal derivatives which would lead to ghosts are not considered. This leads unavoidably to the 
breaking of Lorentz inv ariance at short dist ances. Static and spherically symmetric black hole solutions have been 
studied in this theory KSOOl lLMP09l |P09l |. Considering the Hofava-Lifshitz theory as a modification of General 
Relativity, one can study the solutions of the geodesic equation in the Hoi^ava-Lifshitz black hole space-times. In this 
paper, we are mainly interested in one of the black hole space-times given in [LMPOil] . 

The mathematical techniques described in this paper can not only be used to solve analytically the geodesic equation 
in the Hofava-Lifshitz space-time considered here. The differentials of the first and third kind with underlying 
polynomial curves of arbitrary genus appear in the geodesic equations in many general relativistic space-times. The 
holomorphic differentials appear in the equations for the r- and z9-coordinates, while the differentials of the third kind 
appear in the equations for the ip - and ^-coordinates. This is also the case e.g. in the geodesic equations for neutral 
particles in Taub-N UT lKKHL10j | space-times and in the space-times of Schwarzschild an d Kerr black holes pierced 
by a cosmic string HHLSloj . as well as for charged particles in the Reissner- Nordstrom [GK1C| spa ce-time , where 
elliptic integrals of the first an d third kind appear. They also appear in the Schwarzschild-de Sitter HLOSal . [HLOS | 
and Kerr-d e Sitter sp ace-times HKKL09a( as well as in generalized black hole Plebanski-Demiaiiski space-times in 4 
dimensions HKKL09j with underlying hyperelliptic curves of genus two in the geodesic equations. Also in the higher 
dimen sional space-times of S chwarzschild, Schwarzschild-de Sitter, Reissner-Nordstrom and Reissner-Nordstrom-de 
Sitter HKKLOSi IEHKKLTi | this powerful mathematics of the theory of hyperelliptic functions of higher genera is 
successfully applicable. Geodesies in higher dimensional axially symmetric spac e-tim es, t he Myers-Pe rry space-times, 
are integrated by the hyperelliptic functions of arbitrary genus as well [EHKKLTlj . In |EHKKL11| the integration 
of holomorphic integrals for any genus of the underlying hyperelliptic polynomial curve has been presented. Here we 
expand our considerations and present the solution for the integrals of the third kind for arbitrary genus as well. 



C. Outline of the paper 



The paper is organized as follows. Section |TT] represents a short introduction to the theory of hyperelliptic functions 
that is adjusted to the aim of the paper. In this section we develop the cr-functional realizations of hyperelliptic 
functions. The key- formula in Lemma III. 51 relates the integral of the second kind to the C- function and the 3- vector. 
In Section IIIII we consider the inversion of the holomorphic integral, integrals of the second kind and third kind, 
and also their arbitrary combination. Section IIVI shows that already developed means of computer algebra, like 
Maple/algcurves, are sufficient to compute all period matrices relevant to the cr-functional approach. As a particular 
feature the explicit knowledge of the homology basis ciphered in the software is not necessary. Sections IVllVIII 
exemplify the developed method in the case of genus two, three and four, correspondingly. Section [Villi is devoted to 
the application of the developed method to the problem of geodesic motion in Hofava-Lifshitz black hole space-times. 
We conclude in Section HXl 



II. HYPERELLIPTIC ABELIAN FUNCTIONS 



A. Abelian differentials and their periods 



We first introduce a canonical homology basis of cycles (ai, . . . , a^; bi, . . . , bg), n aj — bi n bj = 0, n bj ~ 
— bi n aj — Sij, where Sij is the Kronecker symbol and n denotes the intersection of cycles. We denote by du{P) — 
(dui(P), . . . , dug{P))^ the basis set of holomorphic differentials (of the first kind) 

dui = dx, j = I,...,5, (II. I) 
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and by dr{P) = (dri(P), . . . ,drg{P))^ the associated meromorphic differentials (of the second kind) with a unique 
pole at infinity, 

dri= ^ (k + 1 - i)Xk+i+i—dx, i = l,...,g, (II. 2) 

k=i ^ 

where the coefficients are as in Eq. p.ip . These canonical holomorphic differentials du and associated meromorphic 
differentials dr are chosen in such a way that their g x g period matrices in the fixed homology basis 

2uj = ( (p dui j , 2uj' = ( <i> du, 

^Ja^ ^ hk=l,...,g VJbfc /i,fc=l,...,S ^^^^^ 



2r/ =(- i dr,) , 2ri' = (- i dr,) 



satisfy the generalized Legendre relation 



with 



MJM^ = -yj (II.4) 



where Og and Ig are the zero and unit g x g-matrices. 

The differential of the third kind flp-^ p^{P) with poles at finite points Pi = (oi, j/i) and P2 = (02, j/2) and residues 
+1 and —1, respectively, can be given in the forrrQ 

o fu\ y + yi dx y + y2 dx 
2(x - ai) y 2[x- a2) y 

If the poles Pi,P2 have the same x-coordinate and lie on different sheets, i.e. Pi — {a,y{a)) and P2 = (a, — y(a)), 
then (jll.6p takes the form 

np^MP) = ^-- (11.7) 

X ~ a y 



The differentials dufc, drk and Q.p-^,p^{P) given above describe the entries in the relation (11.31) . 

We introduce the fundamental bi-differential fl{Q, S) on Xg x Xg which is uniquely defined by the following 
conditions: 

1. It is symmetric, n{Q, S) = n{S, Q). 

2. It has the poles along the diagonal Q — S, namely, if £,{Q) and £^{S) are local coordinates of the points Q and 
S in the vicinity of the point P (^(P) = 0) then the following expansion is valid 

"(^'^) = 7^?7r^Trl^+ E ^rnniPmy'-'aSr-'dmdaS), (11.8) 
(W) t.[^)) ,,„„>! 

where ri„i„(P) are holomorphic in P. 

3. It is normalized such that 



n{Q,S)^0, j = l,...,g. (II.9) 

From (|ir9)) and the bilinear Riemann relation (see, e.g., |BBEIM94 iGHMTOSl ] for details) it follows that the 
b-periods oi fl{Q, S) are 

(f niQ,S)^2i7rdv,iS), j = l,...,5, (11.10) 

Jb-j 

where dv{S) = (dui(S'), . . . , {dvg{S))'^ = (2a;)^^dM(S') is the vector of normalized holomorphic differentials. 



^ One can add an arbitrary combination of holomorphic differentials. However, we take this form as the most simple one which is sufficient 
for the following derivations. 
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We present here the algebraic construction of the fundamental bi-differential fl{Q, S). To do that we will construct 
at first a non-normalized bi-differential T{Q, S) subject to the first two items in the definition of fl{Q, S) above. 

Lemma II. 1. A symmetric bi-differential with the only second order pole along the diagonal defined up to a bilinear 
symmetric form in holomorphic differentials is given by 

r(^'Q) - |-:#^^+dr(z,u;rd«(x,y), (11.11) 
oz Z(x — z) y 

where P ~ {x,y) and Q ~ {z,w), or, equivalently, 

Here dr is the vector of meromorphic differentials (|II.2|) and F{x, z) is a so-called Kleinian 2-polar given by 

a 

F{x, z)^Y. (2^2fe + \2k+i (z + x)) (11.13) 

such that F{x, x) — 2y^ and F{x, z) — F{z,x). 

Proof. Consider the differential of the third kind (jll.6|) 

O fp,_y±w^dx y + w' dx 
2{x — z) y 2[x — z') y 

depending on the variable P — {x, y) and possessing poles in the points Q — (z, w) and Q' = (z', w'). The bi-differential 

tu\A d y + w dxdz m 

-f^g.g,(P)dz = -^^-— (11.15) 

as a form in P has a second order pole along the diagonal P ~ Q, but as a form in Q it has unwanted poles at 
z = oo. This 2-form can be symmetrized (r(P, Q) = T{Q,P)) by adding an additional term dr{z,w)'^du{x,y) that 
annihilates the aforementioned poles. □ 

The differential (jll.lip is defined up to a holomorphic 2-form 2dM"'"(z, w)>cdu{x, y), where x is a symmetric g x g- 
matrix, x'^ = h. This fact will be used for the symmetrization and normalization of the bi-differential. Thus, the 
bi-differential r(P, Q) turns into 

f7(P,Q) = — -f -^dr{z,wfdu{x,y)+2du^{z,w)xduix,y) (11.16) 

oz 2(x — z) y 
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or equivalently 

mQ) - ""^If + 'r -- + ^^-^M>^M^,y) ■ (11.17) 

4(x — zy y w 

The matrix x is chosen so that it norniahzes ^{P, Q) according to (III. 91) and pi.lOp and the factor 2 in the second 
term of pi.l7p is chosen to get precisely the Weierstrafi definitions in the case g — 1. This defines the matrix x in 
terms of the 277- and 2cj-periods as 

X = r/(2w)-i . (11.18) 

We denote by Jac(Xg) the Jacobian of the curve Xg, i.e., the factor C^/F, where F = 2aj ® 2uj' is the lattice 
generated by the periods of the canonical holomorphic differentials. Any point u £ Jac(Xg) can be represented in the 
form 

u = 2uje + 2uj'e' , (11.19) 
where e, e' e R^. The vectors e and e' combine to a 2 x g matrix and form the characteristic e of the point u, 

If It is a half-period, then all entries of the characteristic e are equal to ^ or 0. 

Beside the canonic holomorphic differentials du we will also consider the normalized holomorphic differentials 
defined by 



dv = {2uj)-^du. (11.21) 



Their corresponding holomorphic periods are Ig and r, where the Riemann period matrix t :— ui ^uj' is in the Siegel 



upper half space 6g of g x (^-matrices (or half space of degree g), 

6g = {r g X g matrix|T"^ = r, Im(T) positive definite} . (n.22) 

The corresponding Jacobian is introduced as 

Jac(Xg) (2a;)-ijac(Xg) = CVlg r . (n.23) 

We will use both versions: the first one (2ci;, 2w') in the context of the cr-functions, and the second one (lg,T) in the 
case of the 6'-functions. 

The Abel map 21 : (^g)" — > with the base point Pq relates the set of points (Pi, . . . , Pat) which are called the 
divisor I?, with a point in the Jacobian Jac(Xg) 

^(Pi,...,Pn) V / du. (n.24) 

The divisor V in (in.24p can also be denoted as Pi + ... + Pat - A^Pq. 

More generally the divisor T) is the formal sum 2? = rtiPi + . . . + unPn with integers rij, j = 1, . . . ,N and TV e N. 
The degree of the divisor, deg(I?), is the sum deg(2?) — ni + . . . + un. The divisor of a meromorphic function is 
of degree zero. Two divisors T) and T)' are linearly equivalent if I? — I?' is the divisor of a meromorphic function. 
Linearly equivalent divisors constitute a class. In particular, the canonical class /Cxg is the divisor class of abelian 
differentials of degree 2g — 2. 

The divisor is positive if all Uj > 0. Let 1(1)) be the dimension of the space of meromorphic functions that have 
poles in the points of V of multiplicities not higher than the multiplicity of these points in V. If deg(2?) > g for a 
divisor in general position the dimension 1(7)) is given by 

l(V) ^ deg(V) - g + 1 . (n.25) 

Such divisors are ca lled no n-special. All remaining divisors with deg(T>) > 2 are called special. For a detailed 
explanation see e.g. FK8C| . 
Analogously we define 

2l(Pi,...,P„) = V / d-y = (2c^)-i2t(Pi,...,P„). (11.26) 

k=l -^Po 

In the context of our consideration we choose Pq at infinity, Pq — (00, 00). 
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B. 6— and cr— functions 



The hyperelliptic ^-function with characteristic e is a mapping : Jac(Xg) x 6g — C defined through the Fourier 
series 

6'[e](i;|T) := ^ ^■ni{{m+E')'^r{m+e')+2(v+e)'^ (m+e')} ^ (11.27) 

It possesses the periodicity property 

e[e]{v + n + Tn'\T) = e-2i™'^(''+|rr.')g2i7^(r^^e'-r.'^e)g|J^J(^|^) gj_28) 

For vanishing characteristic we abbreviate 9{v) :— 6[Q\{v\t). 

In the following, the values 6^, e'j, are either or i. The property (|II.28I) implies 

e[e\{-v\T) ^ c-^'''^'^^'e[e\{v\T), (11.29) 

so that the function 6'[e](i;|r) with characteristic e of only half-integers is even if 4e"^e' is an even integer, and odd 
otherwise. Correspondingly, e is called even or odd, and among the 4P half-integer characteristics there are ^(4^ + 2^) 
even and ^(4^ — 2^) odd characteristics. 

The nonvanishing values of the ^-functions with half-integer characteristics and their derivatives are called 9- 
constants and are denoted as 



e[e] := ^?[e](0;T), ^^0[e](^;r) 



, etc. for even [e]; 

z=0 



etc. for odd [e] 



Even characteristics e are called nonsingular if 6'[e] ^ 0, and odd characteristics e are called nonsingular if 9i[e\ ^ 
for at least one index i. 

We identify each branch point ej of the curve Xg with a vector 

/■(ej,0) 

21^ :=/ dit=: 2a;£j -|-2w'e^- e Jac(Xg), j = 1, . . . , 2.g + 2 , (11.30) 

CO 

which defines the two vectors £j and e'y Evidently, [2l2g+2] = [0] = 0. 

In terms of the 2g -f 2 characteristics [2ti] all 4^ half integer characteristics e can be constructed as follows. There is 
a one-to -one co rrespon dence be tween these e and the partitions of the set ^ = {1, . . . , 2(7 -t- 2} of indices of the branch 
points ( [Fav73l |. p. 13, jBak897 p. 271). The partitions of interest are 



U Jm — {«!, ■ • ■ , ig+l-2m} U {ji, . . . }, (11.31) 

where ra is any integer between and [^^] • The corresponding characteristic is defined by the vector 

g+l-2'm 

A„, = 21,, + £m + T£;„ , (11.32) 

k=l 

where Koo G Jac(Xg) is the vector of Riemann constants with base point oo, which will always be used in the 
argument of the ^-functions, and which is given as a vector in ia,c{Xg) by 

all odd [Sj] 

(see e.g. |FK8(| . p. 305, for a proof). 

It can be seen that characteristics with even m are even, and with odd m are odd. There are ^{^g^i) different 

partitions with m = 0, {^g^i) different partitions with m — 1, and, in general, {g^f^2m) down to (^^j^^) = 2g -I- 2 

partitions if g is even and m = or (^^q*^^) = 1 partitions if g is odd and m = {g + l)/2. One may check that the 
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total number of even (odd) characteristics is indeed 2^^ ^ ±2^ ^. According to the Riemann theorem for the zeros 
of ^-functions Fav73l | . 0{Am + v) vanishes to order to at d = and in particular, the function 6{Kao + v) vanishes 
to order [^] at t> = 0. 

Let us demonstrate, following |FK80f . p. 303, how the set of characteristics [2tfe] = [2tfe], k — 1, . . . , 2(7 + 2 looks 
like in the homology basis shown in Figure [T] Using the notation = \{5ik, ■ ■ ■ jSgkY and for the fc-th column 
vector of the matrix r, we find 



(e2fe,0) 



2t 



2g+l 



21 



2g+2 



2g+l 



21 



2g-l 



21 



2s 



''-\e2fc-l,0) 

(e2s,0) ^ 

+ 1:0) ''"^ 
7 1 



(e2g+l:0) 
(62)=, 0) 



(621-1,0) 

Continuing in the same manner, we get for arbitrary I < k < g 

k 



2k+2\ 



[2t: 



2fc+lJ 







.. 


~0 


1 





\ 




1 .. 


1 


1 









fc 












.. 


~0 


1 





\ 




1 .. 


1 









[21: 



■2g+l\ 



f2t^ 



2g\ 



[21; 



2s-lJ 



1/0 
2 U 1 



1/0 
2 U 1 



1/0 
2 U 1 





1 1 



1 

1 1 



1 

1 



(11.34) 



(11.35) 



and finally 



1 fl 



(11.36) 



The characteristics with even indices, corresponding to the branch points e2n, n = 1, . . . , 5, are odd (except for [212^+2] 
which is zero); the others are even. Therefore in the basis drawn in Figure [T] we get 



fc=i 



2k ■ 



(11.37) 



The forn iula pi. 371) is in accordance with the classical theory where the vector of Riemann constants is defined as 
(see Fay |Fav73( . Eg. (14)) 



Divisor = A-((7-l)Po, 



(11.38) 



where A is the divisor of degree <? — 1 that is the Riemann divisor. In the case considered Pq = 00 and A = 
62 + 64 + . . . + e2g — 00. The calculation of the divisor of the differential Hfc^iC^^ ~ e2k)dx/y leads to the required 
conclusion 2A = /Cx, where /Cjf, is the canonical class. 

The Kleinian cj-function of the hyperelliptic curve Xg is defined over the Jacobian Jac(Xg) as 



n{u-M) ■.^Ce[K^]{{2uj)-^u:T)e^ , 
where the symmetric g x g matrix k is defined in (|IL18p . Here G 3ac{Xg) and 



■D a 



goo 



dlt, 



(11.39) 



(11.40) 
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where V = Pi + ... + Pg is a divisor in the general position. The constant 

^1/4 



''=\l6^)\ n , (11.41) 



l<i<j<2g+l 



and M defined in pi.Sp contains the set of all moduli 2uj, 2uj' and 2r/, 2ri'. In the following we will use the shorter 
notation a{u; M) — <t{u). Sometimes the cr-function (jIL39|) is called fundamental cr-function. 

The multi-variable CT-function (|II.39[) represents a natural generalization of the Weierstrafi CT-function given by 

^' V2a.^(ei-e2)(ei-e3)(e2-e3) ^2c^^ \ } 

where i?i is the standard (^-function. We note that (|IL39I) differs in the case of genus one from the WeierstraB cr- 
function by an exponential factor that appears when the shift on a half period in the ^-argument is taken into account 
in the 6'-characteristics. 

The fundamental cr-function PL39I) possesses the properties 

• It is an entire function on Jac(Xg), 

• It satisfies the two sets of functional equations 

cr(M + 2ujk + 2uj'k'; M) = ^^i^k+r^' k- {u+uk+u' k' ) , j^j^ 

a{u;{^M^f)^a{u;M), 

where 7 g Sp(25, Z), that is, 7 J7"^ = J, and is the matrix M with interchanged submatrices oj' and 77. The 
first of these equations displays the periodicity property, and the second one the modular property. 

• In the vicinity of the origin the power series of a{u) is of the form 

cr(it) = S-^{u) + higher order terms , (11.44) 

where St^{u) are the Schur- Weierstrafi functions associated to the curve Xg and defined on 9 [ui, . . . , Ug) by 
the WeierstraB gap sequence at the infinite branch point. For g > 1 it is always a Weierstrafi point. The partition 
TT = [-Kg, . . . , TTi) is defined by the Weierstra B gap seq uence w = (wi, . . . , Wg) as follows: tt^ = Wg-i+i + i — g. 
Details of the definition are given in see also jEHKKLTl| . As an example we will present here the first 

few functions ST^iu) 

5 = 1: Si{ui) = ui, (11.45) 

5 = 2: 52,i(wi,U2) = ^U2-Mi, (11.46) 

5 = 3: 5'3,2,i(wi)'"2,W3) = ttUs - ^U2wi - "2 + ""iws, (11.47) 

45 3 

5 = 4: S'4,3,2,i(ui,U2,W3,M4) = 737:7^1° - -^ulua + —M2ul - "4"! - o"!""! (11.48) 

4725 105 15 3 

The partitions constructed by the WeierstraB gap sequences are denoted in the subscripts. In particular, in the case 
of genus 5 = 4 the partition tt = (1, 2, 3, 4) corresponds to the gap sequence 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, . . ., where 
orders of existing functions are overlined. These are so-called non-gap numbers, in contrast to the gap-numbers. The 
genus is defined by the number of gaps or, equivalently, the first number starting from which no gap appears equals 
2g (in this example 8 = 217). 

The Kleinian Q and p-functions are a natural generalization of the Weierstrafi Q and p-functions and are given by 
the logarithmic derivatives of cr, 

d 

Ci{u) = TT- Incr(M), 

Pijk{u) = -^-^ — ^ — Incr(M) , etc., 
ouiOUjOUk 
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where i,j,k G {l,...,g}. In this notation the Weierstrafi p-function is pii(u). For convenience, we introduce the 
vector of ^-functions C{u) — {Ci(u), . . . , (g{u))'^ and also denote the derivatives of the cr~function by 

d 

CTj(m) = ^cr(M), a,j{u) ^ ———a{u), etc. (11.50) 

OUi OUiOUj 



C. Main formula 



We consider now the integration of the differentials of the second and third kind. 

Proposition II. 2. Let Xg be a hyperelliptic curve of genus g with branch point at infinity. Let D — (Pi,...Pg), 
V ^ {P[,...P'g), Pk = iZk,Wk), P'k = {K^Wl?j be non-special divisors of degree g. Let P = {x, y) and P' = (x', y') 
be two arbitrary points of X . Then 

P g rP^ F{x, z) + 2yw dxdz ^ (£ ' d") ^ ^ (/p„ d« - /^^ du) 

p'hJpL yu, \(;-dn-/-d«) d.) ■ 

Proof. We introduce non-special divisors 

V = ((Zl, W-l), . . . , {Zg,Wg)), V - {{Z[,W'{), {Z'g, W'g)) 

as well as two arbitrary points P ~ {x, y), P' — {x' , y'). The integration of Q) given by (|II.16I) 

V / niP,Q) (11.52) 

yields, according to the Riemann vanishing theorem, 

- ELr It""'^ + Kp„) (/;; d. - ELr It""'' d. + K,^) 

(£ d. - ELr i^:-"^^ d. + i.p„) . (/;; d. - ELr i^:-"^^ d. + Kp„) ' 

where Pq is the base point of the Abel map (that we suppose to be infinity) and Kp^ is the vector of Riemann 
constants with the base point Pq . Using the definition of the fundamental cr-function (|II.39|) we get pi.51[) . □ 

The following corollaries follow from the main formula ()II.51| . 

Corollary II. 3. Let P = {x,y), Pk ^ {xk,yk),k = 1, . . . ,g. Then 



p 

Pa 



F{x, Xk) + lyyk 

Proof. Take the partial derivative jdxdxk on both sides of (|II.51|) . □ 



Corollary II. 4. Let V = Pi + . . . + Pg be a positive divisor of degree g. Then the standard Jacobi inversion problem 
given by the equations 



Pk 



J2 du^u (11.55) 



fe=i 

is solved in terms of Kleinian p-functions a^ 

X^ - Pgg{u)x3-^ - Pg,g-l{u)x3-^ - ... - P g .^{u) = , (II.56) 

yfc = Psss(M)a;r^ + p99,9-iNa;r^ + •■• + P99,i("): k = \,...,g. (11.57) 



In the previous work [EHKKLllI in this formula numbered as (3.44) as well in its particular cases (5.5) and (6.5) the sign "-" was 
misplaced. 
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Proof. To prove PL56I) consider x ^ oo. Substitute x = 1/^^ into (jll.54p . Comparison of the coefficients of 1/^^*^^ 
on the RHS and LHS of (|lL54)) gives relation (|lL56]) . Formula (|lL57)) follows from (|lL56]) and (|lL54| . For details 



see 



BEL93]. 



□ 



For t he me romorphic differentials of the second kind the following relation was proved by Buchstaber and 
Leykin iBLOSj : 



Lemma II. 5 (C- formula). Let Vq be a divisor supported by g branch points e^^ 



, . . . , 



Bi such that 



E 

k=l 



du = Koo : [Kc 



(11.58) 



and V is a non-special divisor of degree g, V = Pi + . . . + Pg. Then for any vector 



the following relation is valid 



u = du e Jac(X) 

J gco 

/ dr = -C(M) + 2(r7V + »7£) + -3(M), 



(11.59) 



(11.60) 



where the components 3j(u) of the vector '^{u) are 

3g(M) = 0, 'ig-l{u) ^ Pggg{u) 

and the other components atl<j<g — 1 are given by the j x j determinants 



3,(«) = 



Pggi'^) 
2pg_i,g(M) 



-1 

Pggiu) -1 



{g - k)pk+l.g{u) Pk+2.g{u) 

(.9 -3 - 1)Pj+2,3(m) p3+3,g{u) 
{g - i)Pi+l,g,g{'^) Pi+2,g,g{u) 



P3s(w) -1 



(11.61) 



The corresponding formula in Buchstaber and Leykin |BL05| coincides with those given by Ba ker 'Ba kslj, P-321, 
only for j = g and j = g — ^ and differs for 1 < j < g — 1. Our further consideration is based on (jll.61[) . 

For our purposes it is necessary to rewrite the expressions for 3j in terms of symmetric functions of the divisor points 
Pi = {xi,yi), ■ ■ ■ ,Pg = {xg,yg). With the solution of the Jacobi inversion problem (IIL56P - (jII.57p the components 
3j{u) yield: 



39-0 

3g-l 



yi 



3g-2 = yi 



{Xi - X2)-- - {Xi - Xg) 
Xi- {X2 + Xg) 



[xi - X2)-- - {Xl - Xg) 



permutations, 
+ permutations, 



)g-3 = yi- 



Xl{X2 + ■ ■ ■ + Xg) + X2X3 + . . . + Xg^lXg 

{Xl - X2)- - ■ {Xl - Xg) 



permutations 



Note that the characteristics in the formula (jll.60|) are not reduced. 

From pLM]) one obtains the relation between the periods of holomorphic and meromorphic integrals 

mk = (Qi^k + J^oo)),.fe^i,...,, , vlk = {Qi^'k + Koo)),.^,, (11.62) 



with 2oJk being the fc-th column of the matrix 2a;. 
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Proposition II. 6. Using (jII.17|) the formula (|IL51I) can he rewritten in the form suitable for the inversion of the 
integral of the third kind 



p 9 r 



fe=i 



y + Wk y + Wl 



P /••D a i du - du 

dit {x,y) / dr{z,w)+\n ) — 

■J-D' ai f du~ f du 

V J OO J QOO 



In- 



CT { du — du 



CT ( f ^ du. — f ^ du 



(11.63) 



where the expression J^, dr{z,w) can be calculated by (|II.60[) . 

D. Stratification of the ^-divisor and inversion 

The 6'-divisor 6 is defined as the subset of Jac(Xg) that nullifies the 6'-function and, therefore, the CT-function, i.e. 

e = {i; e J^c{Xg) I 9{v) = o} . (11.64) 
The subset 6^ C 6, < fc < is called fc-th stratum if each point v E Q admits a parametrization 



dv + Ka 



(11.65) 



where Qq = {Koo} and 6g-i = &■ We furthermore denote 9g = Jac(Xg) and we have the natural embedding 

Go C 91 C . . . C eg_l C Gg = J^(Xg) . (11.66) 

We define the 0-function to be vanishing to the order m(Qk) along the stratum Qk if for all sets a^, j = 1, . . . ,.g 
with < ai + . 



Ug <m holds 



Qai + ...+ag 

du-^' ...du. 



(11.67) 



and there is a certain set of aj, with ai + . . . + ag — m such that (|II.67[) does not hold. The orders m(9fc) of the 
vanishing of 9{Qk + v) along the stratum 9fe for some genera are given in Table H] 

In the following we focus on the stratum 9i corresponding to the variety Qi C Jac(Xg), which is the image of the 
curve inside the Jacobian, 



dv + Ka 



(11.68) 



We remark that another stratification was introduced in |Vanh95| for hyperelliptic curves of even order with two 
infinite points oo+ and oo_ that was implemented for studying the poles of functions on Jacobians of these curves. 
The same problem relevant to strata of the 0-divisor was studied in [AFOOj . 



III. INVERSION OF HYPERELLIPTIC INTEGRALS OF HIGHER GENERA 



A. Inversion of holomorphic integrals 



For the case of genus two the inversion of a holomorphic hyperellipti c integ ral by the method of restriction to the 
-divisor was obtained independently by Grant |Gra90j | and Jorgenson |Jor92j in the form 



(III.1) 



(t(u)=0 
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9 


m(eo) 


m(ei) 


m(e2) 


m(e3) 


771(64) 


777(65) 


771(66) 


1 


1 





- 


- 


- 


- 


- 


2 


1 


1 





- 


- 


- 


- 


3 


2 


1 


1 











4 


2 


2 


1 


1 









5 


3 


2 


2 


1 


1 







6 


3 


3 


2 


2 


1 


1 






TABLE I. Orders 7n(6fc) of zeros 9{Qh + v) at u = on the strata 6^. 



This result was implemented in [E PR03j |. and explicitly worked out in the series of publications |HL08l . iHLOSal . 
| HKKL08 , HKKL09 . HKKLOQ aJ , and others. 

The case of genus three was studied by Onishi 6ni98| . where the inversion formula is given in the form 

CTl3(w) 



0-23 (w) 



M=(U1,M2,U3) • (III.2) 
(t{u)— o-3{u)— 



Formula (jIII.2[) is based on the detailed analysis of the genus three KdV hierarchy and its restriction to the 6'-divisor. 
Below we will present the generalization of pil.ip and pil.2p to higher genera. For doing this we first analyze the 
Schur-Weierstraf5 polynomials that represent the first term of the expansion of (j{u) in the vicinity of the origin it ^ 0. 
The 0-divisor Q and its strata Qk in the vicinity of the origin it ~ are given as polynomials in u. 
An analysis of the Schur-Weierstrafi polynomials leads to 

Proposition III.l. The following statements are valid for the Schur-Weierstrafi polynomials St^{u) associated with 
a partition tt: 

1. In the vicinity of the origin, an element u of the first stratum Qi d Q is singled out by 

SM^O, 7riSn{u)^0 Vj = l,...,5-2. (III.3) 

OUg 



2. The derivatives fulfil 



3. The following equalities are valid for u Cz Qi 



qM QM+1 



(III.5) 



where M = \{g - 2){g - 3) + 1. 
J^. The order of vanishing of restricted to Oi is the rank of the partition tt. 

It was noted in BEL99j | that the Schur-Weierstraf5 polynomials respect all statements of the Riemann singularity 
theorem. In particular, if 



2g-l z^fe-l 



(III.6) 



V25 - I''"' 2fc - 1 ■"' 3 

and if tt is the partition at the infinite Weierstrafi point of the hyperelliptic curve Xg of genus g, then the function 

G{z) -.^ S^{Z - u) (III.7) 



14 



either has g zeros or vanishes identically. Moreover we will conjecture here that the properties of the Schur-Weierstrafi 
polynomials given in Proposition IIII. 1 1 can be "lifted" to the fundamental ct- function (|II.39[) . 

The above ai ialysis perrn its to conjecture the following inversion formula for the general case of hyperelliptic curves 
of genus g>2 [EHKKL11| 



qm+i 

duidu^'^ 



u£0i 



and 



M 



(3-2)(.g-3) 



(Ill.i 



61 = <^ M e Jac(Xg) a{u) = 0, ra{u) ^0 Vj = 1 

I dul 



(IIL9) 



The analog of this formula for strat a Qk, 1 < k < q and (n, s)-curves in the terminology of BEL99j | was recently 
considered by Matsutani and Previato jMP08j |. jMPlo| . 



B. Inversion of the integrals of the second kind 



Also the meromorphic integrals can be expressed in terms of the cr-functions restricted to the stratum 0i of the 
6'-divisor. To do that we consider the ^-formula pi.GOp and move to infinity the g —1 points of the divisor T) on both 
sides of the equality. The poles on both sides cancel and the non- vanishing term s yield the required formula. We 
demonstrate this procedure here for the example of the genus three case given in jEGQ3j | (the analysis of the genus 
two case can be found in EPR03]). For that we use the following proposition 



Proposition III. 2. Let he the genus three curve = 4a;^ + \^x^ + . . 
Then the following formulae are valid 



-Aq. Let u e Oi — {u\a{u) — (73(14) = 0}. 



drs 



dro 



C3 



d?"! = - 



g23(w) _^ 
(72 (m) 

1 q-22(M) 

2 CT2(m) 

1 ai{u)a22{u 



C2 



(III.IO) 



0-12 (m) 
(72 {u) 



Ci 



where P = (x, y) is an arbitrary point of the curve, and Pq — (xq, yo) 7^ (00, 00) is any fixed point. It is convenient to 
choose in particular Pq = (e2g,0). The constants Ci are fixed by the requirement that the right hand side vanishes at 
P = P„ 



Proof. In the case of g = 3 the vector 3(m) (jll.61|) yields 

/3i(«)\ 
3(m) 32(m) = 

V 33(«) / 

First we restrict relation (III.60p to the stratum 02. To do that we use the following expansions 



P33(m)P333(m) + 2p233(M) 
P333(m) 





(III.ll) 







/ dtt 




' 00 


^3 = i/e \ 



If - 



dr 



(£6,0) 




(III.12) 



taking into account the condition <t{u) — in the expansions. Then we restrict the resulting formulae to the 
stratum Qi. To do that we expand the holomorphic and meromorphic integrals du and J^^^ dr in the vicinity 
of P2 — (00,00) taking into account the condition cr^iu) = 0. From that the relations piI.lOl) follow. For the 
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simplification of the obtained relations we used the formulae from (Oni98l| below, which result from the restriction to 
Oi of the KdV hierarchy associated with the genus three curve 



{a333(w)-2a2(«)}|e, =0 





0-233 (m) r~r + '^iw 



CTi33(m) - 



02(m) 



= . 



(III.13) 



□ 



Cases of higher genera cases can be considered in an analogously way though it is not possible to present a general 
formula. 



C. Inversion of the integral of the third kind 



We consider (|II.63|) when all points Pk, PI. have the same coordinate Z'^. = Zk but Wj. — —Wk and choose the 
divisors 



V = {{Z, W), (64, 0), . . . , (e2g, 0)}, V = {{Z, W), (64, 0), . . . , (e2g, 0)} . 



It is evident that 



V 



dr(z, w) 



v 



{zw) 



dr(z, w) 



(e2,0) 



{z,w') 



dr(z, w) 



(III.14) 



(III.15) 



(e2,0) 



From the C^formula pi.60|) we get 



(£2,0) 
(Z~W) 

(£2,0) 



dr(z, w) — — C 
dr{z, w) = C 



(z,w) 

(62,0) 

{Z,W) 

(62,0) 



du + K^]+ 2{'n'e' + r]£) + \'h{Z, W), 



dw + 



2(r,V + 77£)--3(Z,W^), 



where ^giZ, W) = and ior 1 < j < g we have 



w 



(Ill.lf 



k=0 



The Sk{e) are elementary symmetric functions of order k built on 5 — 1 branch points e4,...,e2g: Sq — 1, 5i 
64 + . . . + e2g etc. Then the solution of the inversion problem for the integral of the third kind (jll.63|) takes the form 



W ^—=2/ du^{x,y) C{ du + K^]-2{ri'e' + rje)--5{Z, 

Jp, X- Z y Jp, [ \J{e2,o) J 2 



In 



-In 



W) 



^ [I^ - /g;^^ du - K, 

-(/rdn+/£yd«-K. 



(III.17) 



Note, that (|III.17I) represents a natural generalization of the known antiderivative for elliptic functions 

du 



p{u) - p{v) 

The inversion procedure for the integral of the third kind 



2uC,{v) + In (t(m — w) — In (t(m + w) 



W 



(62,0) 



1 dz 

X " Z w 



= t. 



(III.18) 



(III.19) 
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is described as follows. Here for simplicity the lower bound is fixed at a branch point, say at P' 
In this case the formula (|III.17p takes the form 



W 



62,0) 



1 dx 
X - Z y 



2(2t2 - It) 



C{v + Koo) - 2(r,V + r,£) - i3(Z, W) 



+m — : — — ^ — In- 



a{u + v- Koc) a{Ql2+v- Koo) ' 



where 



iz,w) 



du , m.2 = 



62,0) 



(62,0) 



Au and tt G 0i, u = 



(62,0). 



(III.20) 



D. Solving the inversion problem 

Now we are in a position to describe the inversion procedure for the integral (jl.2p . The formulae given above permit 
to represent equation (|I.2|) in the form 

J-{x, ui, . . . , Ug, constants) = t . (III. 21) 

The function J- depends on x via the g variables wi, . . . , Ug that are abelian images of (a;, y) and various constants 
given by the poles of the integrals of the third kind and the coefficients ak,bk,Ck, and possibly via the elementary 
function £{x) in the equation p.3[) . This relation is complemented by the g — 1 conditions m e Oi 

a(«)-0, —a{u)=0 Vj = l,...,g-2 . (III.22) 

du-'g 

From these relations one can find (numerically) the functions ui = ui{t), . . . ,Ug — Ug{t) and plug these into (|III.8p to 
find X = x{t). 



IV. COMPUTER ALGEBRA SUPPORTING THE METHOD 



Presently effective means of computer algebra are developed to execute the above claimed program of integral 
inversion. We will consider for this, e.g., the Maple/algcurves code. 



A. Riemann period matrix and winding vectors 



For a given curve of genus g we compute first the period matrices (2a;, 2a;') and r = a; ^oj' by means of the 
Maple/algcurves code. From that we determine the winding vectors, i.e. the columns of the inverse matrix, 

(2c.)-i = (J7i,...,J7g) . (IV.l) 



B. Homology basis 

In our analysis we used a specific homology basis for a hyperelliptic curve (see Fig[T]). That was done just to 
clarify the approach. But the result should be independent of the choice of the homology basis. It is possible to 
perform all the calculations without making an explicit plot of the homology basis and to use that one given in the 
Tretkoff-Tretkoff construction that is programmed in the Maple/algcurves. Nevertheless, in certain cases we should 
know the correspondence between the branch points and half-periods in the homology basis used by Maple/algcurves. 
In particular the proposed method of inversion supposes the knowledge of this correspondence. One can find this 
correspondence using the generalized Bolza formulae. That can be done as follows. 

We first find all nonsingular odd characteristics by direct computation of all odd 6'-constants. According to Table U 
we have two sets Bi C 6g-i and B2 C Og_2 of nonsingular odd half-periods. For each element of bi E Bi there are 
Cij^ , . . . , ei^_ J ^ 00 such that 

61 = / dV + ...+ dV + Kooe Sg-l (IV.2) 

7 00 ^00 
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and for each element of 62 € B2 there are e^^ , . . . , ei^_^ 7^ 00 such that 

62 = / dv + ...+ dv + K^e Qg-2 . (IV.3) 

Using the known values of the winding vectors ( |EHKKLll| . Proposition 4.3) one can find the correspondence 
between the sets {ci^, . . . , eijj_^} and {ci^ , . . . , eig_2} of branch points and the nonsingular odd characteristics 
[{2uj)~^ (2lii,...,ig„i) + Kao] and [(2a;)~^ (^ii,...,ig-2) + -K'oo] • Then one can add these characteristics and find 
the one-to-one correspondence 



(%-i.O) 

dvt^ [a,^_J, z = l,...,2g + 2. (IV.4) 



C. Second period matrix 

We present the formula that permits to express the second period matices 21], 2ry' and >r in terms of the first period 
matrices 2uj, 2uj\ the branch points and the ^-constants 

Proposition IV. 1. Let + 2ujKao be an arbitrary even nonsingular half-period corresponding to the g branch 
points of the set of indices Iq = . . . , ig}. We define the symmetric g 'x g matrices 

:= {p^A^x„\,=,_g (IV.5) 

and the g x g matrix H which is expressible in terms of the even non-singular 6-constants^ 

1 



H{^i,) = ^ , where e ^ [{2u:)-^^x, + . 

Then the x-matrix is given by 



The half-periods rj and rf of the meromorphic differentials can be represented as 



ri = 2xio, 77' = 2W-y(a;~^)'^. (IV.8) 
We remark that pV.7l) represents the natural generalization of the Weierstrafi formulae^ 

see e.g. the Weierstrafi-Schwarz lectures, |Wei893l | p. 44. From (HVJl) follows 



(IV. I 



-i^(2lxo) - i((2a.)-i)^i7(2lxo)(2c.)-i . (IV7) 



12 MO) MO) 

Therefore Proposition IIV.II allows the reduction of the variety of moduli necessary for the calculation of the a- and 
p-functions to the first period matrix. The generalization of pV.lOp to arbitrary higher genera algebraic curves has 
recently been discussed in KShll] where its role in the construction of invariant generalizations of the Weierstrafi 
cr-function to higher genera is elucidated. 

For the following Lemma we introduce g vectors 

= (l,e,,...,ef , jelo. 



In the previous work lEHKKLllI this formula numbered as (3.50) was written in an equivalent but more complicated form with misplaced 
sign "-" . 

* The correspondence between the numeration of branch points and S-constants i9j (0) depends on the chosen homology basis 
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Lemma IV. 2. Let Sj! he elementary symmetric functions of order k built on the 2g—l elements {ei, ... ,e2g+i}/{ei,ej}. 
Then the following relations are valid 

^ * ^ ^ n— 

Here F(x,z) is the Kleinian 2-polar (|II.13[) . 

Proof The LHS of the relation pV.lip fohows from (llL54t 

where the argument of the p-function is a non-singular even half period given as an abelian image of g branch points 
Pk = (cifc ,0), ii, . . . ,ig G {1, . . . , 2g + 1}. The branch points and on the right hand side of pV.12[) belong to 
the chosen subset of g branch points. The RHS of the relation pV.ll[) can be checked directly for small genera and 
presents a conjecture for higher genera. 

□ 

Varying the integers i,j along the set Xq one can obtain from (IV. 10) g{g — l)/2 equations with respect g{g — l)/2 
components, pg_i^g_i, . . . , pn of the matrix and then solve them by the Kramer rule. Indeed, since the 

components pgi(Io) are already known from the solution of the Jacobi inversion problem (11.58) then Eq. (IV. 10) can 
be simplified for every pair i ^ j G Tp as follows 

\ keXo/{i.j} J »=0 

where 

$(2liJ = (p.,(2li„)),,,=i,...,g„i, E.,^{l,e,,...,erY ■ (IV.14) 

Examples for genus two, three and four in Sections IVllVIII show that the entries in *p are polynomials in the branch 
points Ci (see Eqns. (jV.14[) . (jVI.19|) and (IVII.21I) '). We infer that this is true in general but do not present here the 
general expression of these polynomials. 

D. Characteristics 

We explained above how to find the correspondence between the half-periods and the abelian images of the branch 
points on the basis of Bolza type formulae. This correspondence is necessary for the description of the real evolution 
of the system that corresponds to the motion of the divisor points, i.e. functions on the upper bound in the Abel 
map, over the Riemann surface between certain branch points or between the corresponding half-periods in the Jacobi 
variety. We intend to reduce a number of cases when complete hyperelliptic integrals are calculated. Finding the 
above correspondence allows to present initial/final points of the evolution in terms of a half period. 

E. Vector of Riemann constants 

The vector of Riemann constants Koo enters the derived inversion formulae and can be computed as follows. It 
was proved that the vector of Riemann constants with the base point fixed in a branch point, e.g. infinity in our 
considerations, is a half-period. From Table |T] for the stratum Qq follows whether K^o is even or odd (parity of m) 
and whether it is singular (m > 1). In accordance with the definition of the 0-divisor it is sufficient to find the even 
or odd half-period K^o which satisfies the condition 

9 r d« + - 0, 

\k=i-'°° J 
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where Pi, ... , Pg-i are g — 1 arbitrary points on the curve Xg. 

Alternatively one can use the correspondence between the branch points and the half periods. Among the 2g + 2 
characteristics (|IV.4|) there should be precisely g odd and g + 2 even characteristics. The sum of all odd characteristics 
gives the vector of Riemann constants with the base point at infinity. This characteristic will then be singular of order 

V. HYPERELLIPTIC CURVE OF GENUS TWO 



We consider a hyperelliptic curve X2 of genus two 

= 4(z - ei){z - e2){z - e3)(z - e4)(z - 65) 
= 4z^ + XiZ^ + Xsz^ + \2Z^ + Aiz + Ao . 

From pi.ll) and (jll.2[) the basic holomorphic and meromorphic differentials are 



zdz , 



Then the Jacobi inversion problem for the equations 



is solved according to pi.56p and pi.57p in the form 

Zl+ Z2 = p22{u 

Wk ^ p222{u)zk + pl22iu), fc=l,2. 





/•(Z2,W2 


dz 




'J 00 


w 


/■("^■"'i) zdz 


r(z2.W2] 


zdz 


Joe W 


•/ 00 


w 


in the form 






- Z2 = P22(m), 


Z1Z2 = - 


-Pl2 



Ml , 



U2 



(V.l) 



d.i = ^, d.i = l^^l±^^^^dz, (V.2) 

w Aw 



du2 — , dr2 = — dz . (V.3) 

w w 



(V.4) 



(V.5) 



A. Characteristics in genus two 



The homology basis of the curve is fixed by defining the set of half-periods corresponding to the branch points. The 
characteristics of the abelian images of the branch points are defined as 



[2t. 



(e.,0) 



du 



which can be also written as 

^, = 2uj£, + 2uj'e'i, i = l,...,6. 
In the homology basis given in Figure [2] the characteristics of the branch points are 



[2ti 



1/1 0^ 
2 lo 0, 



1/1 0^ 

Hi Oy 
1/00^ 
2 1 1, 



[2t3 



1 V 

2 U 0, 



The characteristics of the vector of Riemann constants Kao yield 



1 1 1' 

2 1, 



(V.6) 



(V.7) 
(V.8) 

(V.9) 
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ea — oo 



FIG. 2. Homology basis on the Riemann surface of the curve X2 with real branch points ei < 62 < • . • < ee = 00 (upper sheet). 
The cuts are drawn from 62^-1 to 62,, i = 1, 2,3. The b-cycles are completed on the lower sheet (dotted lines). 



From the above characteristics we build 16 half-periods. Denote 10 half-periods for i ^ j = 1, . . . , 5 that are images 
of two branch points as 



= 2oj{£, + £j) + 2uj'{£[ + e;.), i - 1, . . . , 5 . 
Then the characteristics of the 6 half-periods 

[(2t^)-i2l, + Xoo] -:<5n 1 = 1,..., 6 
are nonsingular and odd, whereas the characteristics of the 10 half-periods 

[{2u)-^n,J^'Koo]^■■e^J, l<i<j<5 



(V.IO) 
(V.ll) 
(V.12) 



are nonsingular and even. 

Odd characteristics correspond to partitions {6} U {1, . . . , 5} and {fc} U {ii, . . . , 14, 6} for ii, . . . , 14 7^ k. The first 
partition from these two corresponds to 0o and the second to 61. 

From the solution of the Jacobi inversion problem we obtain for any i, j — 1, . . . , 5, i ^ j 



From the relation 



one can also find 



ei + ej = p22(J^'y), -eiCj = pi2(rJij) . 

F{xi,X2) - 2yiy2 



(V.13) 



Pii(m) 



ej{ep + eq + e,,) -I- epEqCr = pii{ftij) , 



(V.14) 



where ij,p,q, and r are mutually different. 

From (|V.13|) and (jV.14[) we obtain an expression for the matrix x (jIV.7|) that is useful for numeric calculations 
because it reduces the second period matrix to an expression in the first period matrix and ^-constants, namely, in 
the case Ci — ei , ej = 62 , 



1 ( 6162(63 + 64 + 65) + 636465 -6162 \ I, ^_iT 1 / 6'ii[e] 6'i2[e] \ ^_l 

" ~ ~H -^^^^ er 62 j " 2^'"^ m [ e^2[e] ^22[e] ) ^'"^ ' 

where the characteristic e in the fixed homology basis reads 

s = [2li] + [2I2] + [K^] = ( I I j . 



(V.15) 
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B. Inversion of a holomorphic integral 

Taking the limit Z2 ^ oo in the Jacobi inversion problem (|V.4p we obtain 



— / ^U2. V.16 



The same limit in the ratio 

P12(m) Z1Z2 



P22{u) Z1+Z2 

leads to the Grant- Jorgenson formula (jllLlI) . In terms of 0- functions this can be given the form 



z = 



(V.17) 



(V.18) 



?((2w)-1m:t)=0 



where here and below dui = ^ij^ derivative along the direction Ui. The "winding vectors" Ui, U2 

are the column vectors of the inverse matrix (2a;) ~^ (|IV.1[) . 
From (jlll.ll) we obtain for all finite branch points 

^» = ^oTT' or equivalently, e^^-- — -7-ry, i = l,...,5. (V.19) 

This formula was mentioned by Bolza |Bol886j (see his Eq. (6)) for the case of a genus two curve with finite branch 
points. 

The ^-formula reads 



Ci(m) + 2ni + = / dri{z,w)+ / dri{z,w) 

'(e2,0) "'(64,0) 



2 Zi — Z2 



-C2(M) + 2n2- / dr2{z,w)+ dr2{z,w) , (V.20) 

•^(62,0) "'(64,0) 

where Uj = Vji^'i + Vji^i- Here the characteristics e- and Si of Koo are not reduced. Choosing (zi, uji) = {Z, W), 

{z2,W2) = (64,0) we get from (|V.20p 

/ i-{z,w) \ 1 W /•(^.'^) 
-Ci / dw + Koo +2ni + -- =/ dri{z,w) , 

y{e2-0) J 2Z-e4 J(g^o) 

-C2 / dtt + Koc +2n2 = / dr2(z,w;) . (V.21) 

V"'(e2,0) / "'(62,0) 

The inversion formula for the integral of the third kind (IIII.17P is written as 

W r = 2 (u- - u'-) r'^' dr + m--;---^-; - ^^.ju'-v-K^) ^^^^^^ 

x~Zy \ ^ "/(ea.o) a{u + v-Koo) a{u' + v-Koo) 



P' 



with 



V 



(ZM) pP pP' 

du, u — du, u' — I du 

(62,0) J 00 



r(Z,W) 

and ueBijit'eOi. The integrals / dr are given by the formula (jV.21[) 

J(e2.0) 



'(62,0) 

In the case when the base point P' is chosen to be a branch point, say (e2, 0) then the final formula takes the form 
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W 



1 dx 



(62,0) x~ Z y 



{u^ - ml) 



^^ a {u-v- Koo) _ ^^ o-(2t2 - V- Kpg) 



(V.23) 



VI. HYPERELLIPTIC CURVE OF GENUS THREE 



We consider also a hyperelliptic curve Xg, of genus three with seven real zeros given by 



w 



4(z - ei){z - e2){z - ei){z - ei){z - e5)(z - ee){z - 67) 



= Az"^ + Xqz^ + . . . + Aiz + Ao . 
The complete set of holomorphic and meromorphic differentials with a unique pole at infinity is 



dui = — , 
w 

zdz 

dU2 = , 

w 

z'^dz 

dU3 = , 

w 

The Jacobi inversion problem for the equations 



dz 

dn = z{2Qz'^ + AXqz^ + 3A50^ + 2XiZ + A3) — 

Aw 

dz 

dra = ^^(12^2 ^ 2\qz + A5) — , 

Aw 

drs = 



w 



— = Wl, 

w 



U2, 



is solved by 



f^'''""'^ zdz r^'''""'^ zdz ^ /■(^^■"'^) £d£ 
Joo W W W 

/ + / + / = U3 

Joo W Joo W W 



Zi + Z2 + Z3 = p33{u), Z1Z2 + Z1Z3 + Z2Z3 = -p23iu), Z1Z2Z3 ^ pi3{u) 
Wk = p3ii{u)zl + p2is{u)Zk + P133(m), = 1, 2, 3 . 



(VI.1) 



(VI.2) 



(VI.3) 



(VI.4) 



A. Characteristics in genus three 
Let 2lfc be the abelian image of the fc-th branch point, namely 



2t 



k 



du = 2wefe + 2w e^, 



1 



(VI.5) 



where and are column vectors whose entries Skj, e'l^j are ^ or for all A; = 1, ... , 8, j = 1, 2, 3. 

The correspondence between the branch points and the characteristics in the fixed homology basis is given as 



[2li] = 



[2l4 



1 /l 

2 i^O Oy 

1 /o 1 0^ 

2 I 1 1 0; 



[9t2] 



1/100' 
2 1^1 0^ 

1/001^ 
2 1,11 0, 

1/000^ 
2 111, 



[2t3] 



1/010' 
2 1^1 0^ 

1/001^ 

1/000^ 
2 io 0, 



(VI.6) 
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FIG. 3. Homology basis on the Riemann surface of the curve X3 with real branch points ei < 62 < . . . < es = 00 (upper sheet). 
The cuts are drawn from e2i-i to 62;, i = 1, 2,4. The b-cycles are completed on the lower sheet (dotted lines). 



The vector of Riemann constants K^o with the base point at infinity is given in the above basis by the even singular 
characteristics, 



K^] = [2I2] + [2t4] + [2l6] = - r ^ M . (VI.7) 




From the above characteristics the 64 half-periods can be built as follows. If we start with singular even character- 
istics, then there should be only one such characteristic that corresponds to the vector of Riemann constants ifoo- 
The corresponding partition reads J2 = {} U {1, 2, . . . , 8} and the ^-function 6{Koo + v) vanishes at the origin 
V = to the order m = 2. 

The half-periods Ai = (2a;)~^2lfe -|- K^o € Oi correspond to partitions 

XiUJi = {fc,8}U{ji,...,j6}, ji,...,j6^{8,fc} (VI.8) 

and the ^-function ^(Ai -|- v) vanishes at the origin t; = to the order m = 1. 
We also denote the 21 half-periods that are images of two branch points 

flij = 2uj{ei + ej) + 2j{e\ + e'^), i,j = l,...,7,ij^j. (VI.9) 
The half-periods Ai = {2u!)~^^lij + € Q2 correspond to the partitions 

XlU Jl ={z,j}U{jl,...,.j6}, Ji,...,j6 ^ (VI.IO) 

and the ^-function 6{^i-\-v) vanishes at the origin, v = 0, as before to the order m = 1. Therefore the characteristics 
of the 7 half-periods 

[{2uJ)-^^2i, + K^]=■.5,, z = l,...,7 (VI.U) 

are nonsingular and odd as well as the characteristics of the 21 half-periods 

[{2io)-^n,j + Koo] =:e,j, l<i<j<7. (VI.12) 

We finally introduce the 35 half-periods that are images of three branch points 

nijk = 2uj{£i + £j+ek) + 2uj'{ei + £j + £'k)e3a.c{Xg), l<i<j <k<7. (VI.13) 

The half-periods A2 = {2oj)~^Clijk + K^o correspond to the partitions 

XoU Jo = {«,i,fc,8}u{ii,...,i4}, h,---,ki{i,j,k,^}- (VI.14) 

The ^-function 6{A.2 + v) does not vanish at the origin t; = 0. 
Furthermore, the 35 characteristics 

Eijk = [{2uj)-^nijk + ifoo] , l<i<j<k<7 (VI.15) 

are even and nonsingular while the characteristic [Kao] is even and singular. Altogether we got all 64 = 4^ charac- 
teristics classified by the partitions of the branch points. 
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B. Inversion of a holomorphic integral 



All three holomorphic integrals, 



= Ml, 

W 



f^'"') zdz /•(^■"') z^dz 

J oo 



W 



W 



U3 



(VI.16) 



are inverted by the same formula pil.2[) . Nevertheless, there are three different cases for which one of the variables 
Mi,W2,U3 is considered as independent while the remaining two result from solving the divisor conditions (j{u) = 

CT3(m) =0. 

Formula piL2p can be rewritten in terms of 6'-functions as 



where 63 = (0, 0, 1)"^. This represents the solution of the inversion problem. 

From the solution of the Jacobi inversion problem follows for any l<i<j<k<7, 

+ Cj + Cfc = p33(r2ijfe), -e^e-,- - e^efe - ejCk = p23{^ijk), eiCjCk = piai^ijk) ■ 
Solving equations (jIV.lip we find 

Pl2i^ijk) = -8381 - 5*4 , pll{flijk) = ^3^2 + 5154 , p22(i^ijfe) = 5*3 + 2s3 + S2<5'l , 



(VI.17) 



(VI.18) 
(VI.19) 



where the s; are the elementary symmetric functions of order I of the branch points e^, Cj, and Si are the elementary 
symmetric functions of order I of the remaining branch points {1, . . . , 7} \ {i, j, k}. 

From (jIV.7|) using (jVI.lSp and (jVL19[) one can find the expression for the matrix x. To do that we take the 
half-period Jli23 



with 



(VI.20) 



0ii[e] 0i2[e] 013 [e] 
^^(^^123) = -zh I ^12 [e] 022 [e] 023 [e] I ^ £ = 
013 [e] 023 [e] 033 [e] 



i i 

2^2 

i i 
2^2 



(VI.21) 



For the branch points ei, . . . , eg the expression 



is valid. Furthermore we have for i, j — 1, . . . ,8, i ^ j 



and for i = 1 , . . . , 7 



dui 


du3 


+ 2mjxe3 


0[Koo]ii2oo)-^mf,T) 






+ 2Qi[Ke3 


0[Koo]((2w)-121,;t) 


1,... 


,8, I 






ei ^ 


-ej = 


(^3{^i]) 


_ 9i720[£ij] 


BiCj - 


_ ai{n.ij) _ 
a3{flij) 








f7l(2t.) 




e 


i — 


(72(21.) 


^uAS^] ■ 



(VI.22) 



(VI.23) 



(VI.24) 



The C-formula reads 
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1 Wl(zi — Z2 — Z3) 

Ci(u) + 2ni + -7 77 7 + permutations 

2 (21 - Z2)(Zi - Z3) 

(zi.Wl) /•(Z2,W2) p(z3,W3) 

dri{z,w) + / dri{z,w) + / dri(z, i/j) 

(62,0) -/(ei^O) "'(£6,0) 

^ , , „ 1 Wi 

- C2[u) + 2x12 + r-. r + pcrmutations (VI.25) 

2 (Zi - Z2){zi - Z3) ^ ' 

/•(zi.Wl) /•(Z2,W2) p(z3:W3) 

= / dr2{z,w)+ / dr2{z,w)+ / dr2{z,w) 

J (€2,0) "'(64,0) "'(66,0) 

pizi-Wi) p{z2,W2) piz3,W3) 

-C3(M) + 2n3= / dr3{z,w)+ / dr3{z,w) + / dr3{z,w) , 

"'(62,0) "'(64,0) "'(66,0) 

where rij = J^^^iVji^'i + Vji^i- Here the characteristics and e,; of /too are not reduced. Choosing (zi,wi) = 
iZ,W),{z2,W2) = (64,0), (Z3,«^3) = (66,0) wc get from dVL25]) 

- f d„ + ic J + 2,., + ij|<^£|±^ ^ , 

V-/(62,o) / 2(Z-e4)(Z-e6) ^(^^^0) 

-C2 / du + K^]+2n2 + -— r— -=/ dr2{z,w), (VI.26) 

\"'(62,o) / 2(Z-e4)(Z-e6) 7(^2,0) 

-Cs / dit + Koo +2n3= / dr3(z,w) . 

\"'(62,0) / "'(62,0) 

The inversion formula for the integral of the third kind (IIII.17P is written as 

Wr^'^^-2(u-- u'-) r""' dr + In-j---^-; - ^^a(y-v_K^ ^^^^^^^ 
Ip, x-Z y \ / .h„^n\ (7 (u + V - K^,) a iu' + V - Kr^) 



with 



V 



'Z.W) pP pP' 

dit, u = dit, u' — du 



'(62,0) 

f{Z,W) 

and Me6i,M'e6i. The integrals / dr are given by the formula (|VI.27p . 

"'(62,0) 



VII. HYPERELLIPTIC CURVE OF GENUS FOUR 

As the next example we consider the hyperelliptic curve X4 of genus four with nine real zeros given by 

9 

4 [|(z - 6fc) = + Xsz^ + . . . + Aiz + Ao . (VII.l) 

k=l 



All calculations in this section will be done without explicit plotting of the homology basis. 

The complete set of holomorphic and meromorphic differentials with a unique pole at infinity is 

dui = — , dri z(A3 + 2A4Z + SX^z'^ + AXqz^ + 5X7Z^ + 6Xsz^ + 28z^)^ , 

W 4:W 

du2 = — , dr2 = z^(A5 + 2A6Z + SAyz^ + 4A82^ + 20z'^)^ , 

du3 = ^, dr2 =z3(A7 + 2A8Z + 12z2)— , (VII.2) 
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z^dz , z'^dz 
du4 = , drs = 

w w 

The Jacobi inversion problem for the equations 



— = Ui, 

w 



r^'"""'^ zdz /•("^■"'^) £d£ r^^'^^^'' zdz [■^''■'^'^ Zdz _ 

Joo W W W W 

ioo W W 



= Ui 



where Sfe and are column vectors whose entries ek,j, £^ ^ are 5 or for all = 1, ... , 8, j = 1, 2, 3. 
The characteristics of the branch points in the fixed homology basis yield 

. , _ 1 A 0\ ro. 1 - 1 A 0^ fOl 1 - 1 1^0 1 0\ 
^"^^^-2 1^0 oj ' ^"^^^-2 1^1 oj ' ^"^^^-2 1^1 oj ' 



2 1 oy ' ' 2 1 oy ' ' 2 1 1 0^ 

1 /O A ro. 1 - 1 1^ f2l 1 - ^ 0^ 
[^^1-2 1^1 1 1 oj ' [^^1-2 1^1 1 1 ij ' t^^]-2lllll. 

[2lio: 



1 




























1 







1 





















1 


1 





1 /0 

2 io 0, 



(VII.3) 



is solved by 

4 

Y^Zi = P44(m), Z1Z2 + Z1Z3 + Z2Z3 + Z-1Z4 + Z2Z4 + Z3Z4 = -p34{u), 
i=l 

Z1Z2Z3 + Z4Z1Z2 + Z4Z3Z1 + Z4Z3Z2 = p24(w), ^122^3^4 = -Pi4(m) 
Wk = p4Ai{u)zl + p34i{u)zl + p244{u)Zk + pl44(w), fc = 1, . . . , 4 . 

A. Characteristics in genus four 
Let 2tfe be the abelian image of the fc-th branch point, namely 

2lfc=/ du = 2uiek + 2u)'£'i^, A: = l,...,10, (VII.4) 



(VII.5) 



The characteristics of the vector of Riemann constants K^o with the base point at infinity are even and singular 
as in the genus 3 example 

[Koo] = [2I2] + [2I4] + [2i6] + [stg] = ^ I J j) • (vn.6) 

From the above characteristics 256 half-periods can be built as follows. If wc start with singular even characteristics, 
then there should be only one such characteristic that corresponds to the vector of Riemann constants K^o- The 
corresponding partition reads I2IJJ2 = {} U {1, 2, . . . , 10} and the ^-function 6{Koo + v) vanishes at the origin w = 
to the order m = 2. ^ 

The half-periods Ai = (2w)~^2lfe -|- K^o G 61 correspond to partitions 
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XiU Ji = {fc,10}U{ji,...,j8}, ji,...,j8^{10,fc} (VII.7) 

and the ^-function 9{A.i + v) vanishes at the origin v = to the order m = 2 as follows from the Table HI The 
characteristics of the 9 half-periods 

[{2Lj)-'<2i, + K^]=:5,, ^ = l,...,9 (VII.8) 

are singular and even. 

Also denote the 36 half-periods that are images of two branch points 

n,, =2uj{e, + e,) + 2uj'ie'^+e'j), i,j = l,...,9,i^j. (VII.9) 

The half-periods A2 = (2^)^^$!^^ + Kco G ©2 correspond to the partitions 

IiiiJi^{i,j}U{ji,...,js}, Ji,---,j8i{i,j} (VII.IO) 

and the 0-function 6{A.2 + v) vanishes at the origin, v — 0, to the order m = 1. Therefore the characteristics of the 
36 half-periods are nonsingular and odd 

[(2w)-ij^y + K^] =: , 1 < z < j < 9 . (VII.ll) 

We introduce 84 half-periods that are images of three branch points 

riijk = 2uj{ei + Ej + £fe) + 2w'(£- + e'j + e'^) e Jac(Xg), 1 < i < j < fc < 9 , (VII.12) 

and 126 half-periods that are images of four branch points 

n.jki = 2a;(£, + £j +ek + £;) + 2uj'{e', + e'^ + £;,£;) G Jac(Xg), 1 < z < j < /c < 9 . (VII.13) 

The half-periods A3 — (2cL>)^^Jlijfe -I- K^a correspond to the partitions 

lo U Jo = {z, J, k, 10} U {ji, . . . , je}, ji, . . . , is ^ {hJ, fc, 10} . (VII.14) 

The 6*- function 6{^j, + v) vanishes at the origin 7; = to the order m = 1. The 84 characteristics 

e,jk=[{'2oj)-^n,jk + Koo], l<i<j<k<9 (VII.15) 

are odd and nonsingular as follows from the Table HI The half-periods A4 — {2uj)^^Qijki + Koo correspond to the 
partitions 

lo U Jo - {i, J, k, 1, 10} U {ji, . . . , J5}, Ji,---J5i {i, J, k, 1, 10} . (VII.16) 
The 6'-function 6'(A4 + v) does not vanish at the origin v = Q. And the 126 characteristics 

e^jki = [i2u)-'n,jki + K^] , l<i<j <k<9 (VII.17) 

are even and nonsingular. 

Altogether there are 256 = 4'^ characteristics classified by the partitions of the branch points. 

B. Inversion of a holomorphic integral 

For the case of genus four the formula (|III.8P reduces to 

(7144 (m) 



(7244 (m) 

The four holomorphic integrals, 



u = (m1,W2,M3,U4) • (VII.18) 

(t{u) — 0,(74 (u)—0,a44 (it) — 0, 



dz /■(^■"') zdz z'^dz /■(^''") z^dz 

Ul, / =U2, / =U3 / =U4 (VII.19) 



00 



w 1^ w 
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are inverted by the formula (|VIL18p . 

From the solution of the Jacobi inversion problem (|II.56p for any 4 roots l<i<j<k<9 follows 

Ci + Cj + Cfe + e; = p44{^ijki), eiBj + CiCk + ejCfe + e^e; + e-,e/ + efeg; = -p^ii^ijki), 

eiCjek + eieiBj + eietei + e/efcCj = p24(J^ijfe/), eiejCkCi = -piA[^ijki) ■ (VII.20) 
From equation (jIV.ll[) one finds the remaining components of the function pij [u) 

Plli^ijkl) = S2S5 + S4S3 , pui^ijkl) = -546*2 - 515*5 , pisi^ijkl) = 5*5 + S4S1 , 
P22i^ijkl) = 255 + SiSi + S3S2 + 25iS4 , 

p23{^i]kl) = -S35i - 54 - 2S4 , pssi^ijkl) ^83 + S2S1 + 2s3 , (VII. 21) 

where si are the elementary symmetric functions of order I of the branch points e^, Cj, Cfc, e; and Si are the elementary 
symmetric functions of order I of the remaining branch points {1, . . . , 9} \ {«, j, k, I}. 

From pV.7|) using (jVII.20|) and (IVIL2ip one can find the expression for the matrix x. To do that consider the 
half-period Jli234 



-^*P(ni234) - ^(2c.)-1^77(ni234)(2w)-l 



(VII.22) 



with 



1234) 



Mk,=i 



1 1 1 1 

2 2 2 2 

^ i 



For the branch points ei, . . . , eio for g > 3 one can also use the formula jEHKKL"Tl| 



(74 (73 



where T = I" U { i} wit h I" = {ii,..., 12} and i 7^ 10, i ^ I". 
The C-formula (|II.60P for genus four reads 



Ci(M) + 2ni + - 



1 Wi{zl - Zi(z2 + 2:3 + Z4) + Z2Z4 + Z3Z4 + 2223) 



(21 - 22) (21 - 23) (21 - 24) 

(zi.lDl) AZ2,W2) r{z3,W3) 

dri(2,u;)+ / dri(2,u')+ / dri{z,w) 

(62,0) "'(64,0) "'(66,0) 



permutations 

{z4,,Wi) 



dri(2, w) 



(68,0) 



/■(.\^->„ _L 1 W;i(zi-(z2 + 23 + 24)) 
C2(M)+2n2 + - 



2 (21 - 22)(2i - 23)(2i - 24) 

Ar2{z,w) + / dr2(2,w)+ / dr2(2,ii;) 

(62,0) "'(64,0) "'(66,0) 



permutations 

/•(Z3,W3) 



(24,104) 



dr2(2, w) 



(e8,0) 



C3(M) + 2n3 + i 



Wi 



permutations 

n(23,U'3) 



2 (21 - 22)(2l - 23)(2l - 24) 
r(zi,Wi) I'(Z2,W2) 

dr3(2,w)+ / dr3(2,w)+ / dr3(2,w)- 

'(62,0) "'(64,0) "'(66,0) 

r(zi,Wi) p{z2,W2) /■{Z3,W3) 

C4(M) + 2n4= / dr4(2,u))+ / dr4(2,w)+ / dr4(2,u') 

"'(62,0) "^(64,0) "'(66,0) 



(24,104) 



dr3(2, w) 



(68,0) 

(Z4,«'4) 



dr4(2, w) 



(68,0) 



where rxj = Ylt^i Vji^'i + Vji^i- Here the characteristics e'^ and Ei of K^o are not reduced. 

Choosing (21, wi) = (Z, VF), (22, W2) = (e4,0), (23,103) = (egjO), (24,^4) = (e8,0) we get from (|VIL25p 



(VII.23) 



(VII.24) 



(VII.25) 



Ci / du + K, 

V(e2,0} 



2ni H ^ 7 — 

2 {Z - e4){Z - ee){Z 



6466 + 6468 



-C2 



(Z,W) 



du + K^] + 2n2 + 



(62,0) 



1 WjZ-iej + ee + es)) 

2 {Z ~ e4){Z - ee){Z ~ es 



{Z.W) 

(62,0) 

(Z.W) 

(62,0) 



dri(2, vS) 



dr2{z,w) 
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-Cs / du + K^] +2n3 + -— -— -— -=/ dr3(z,^YII.26) 

( r{Z,W) \ AZM) 

-C4 / dtt + Koo +2n4= / dri{z,w) . 

\J{e2,0) J "'(62,0) 

With (|VII.27|) the inversion formula for the integral of the third kind piI.17[) yields 

. -2 - u'-) r^""' d. + m-^j-"--^-; - m-;-'---^-; (vii.27) 

Ip, x-Z V \ / .h^„n\ a iu + V - Kryr,] a lu' + V - Knc) 



with 



r{Z,W) nP pP' 

V = du, u = du, u' = du 

J (e2 ,0) 00 ^00 

and M e 61, It' e 61. 

VIII. APPLICATION: SOLUTIONS TO THE GEODESIC EQUATION IN HORAVA-LIFSHITZ BLACK 

HOLE SPACE-TIMES 



Now we are applying our developed methods of integration of differentials of the first and third kind to the integration 
of the equations of motion of pointike test particles in Hof ava-Lifshitz space-times. This class of space-times provides 
a quantum gravity space-time model which is power-counting renormalizeble and reduces to General Relativity in the 
infrared limit, i.e. at large distances. However, it faces the problem to violate Lorentz-symmetry at short distances. 
The main reason for that is that the model contains only higher order spatial derivatives in the action, while higher 
order temporal derivatives (which would lead to ghost degrees of freedom) do not appear [H09aLlH09bt . 

A. Equations of motion 

The metric of a spherically symmetric black hole in Hoi^ava-Lifshitz gravity is given by 

ds^ = N^{r)dt^ - f-\r)dr^ - {dd^ + sin^ Bdip'^) . (VIII.l) 
The Lagrangian for a point particle moving in this space-time reads 

where e = for massless particles and e = 1 for massive particles, respectively. 

The constants of motion are the energy E and the angular momentum (direction and absolute value) of the particle. 
We choose = 7r/2 to fix the direction of the angular momentum and have 

E:^N'^, L,:^r'^. (VIII.3) 
dr dr 



Using these constants of motion we get 



with the effective potential 



K=ff(r) =iV' (e + ^ ) . (VIII.6) 
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Static spherically symmetric black hole solutions of this theory have been discussed in (KS09I . ILMP09I . |P09l | . In all 
cases, the metric functions are of the form 



N^{r) = f{r) = 1 + cir^ - y/c2r^ + c^r , (VIII.7) 

where ci, C2 and C3 are constants. Here, we will be interested in the case 

ci=-Aw, C2 = 0, C3 = aV-Aw , (VIII.8) 

where Aw is proportional to the negative cosmological constant and a > 4/3'^/^ is a n arbitrary p arameter |LMP09j |. 
The geodesic equations of point particles in the fields given by th e solutions i n KSOOl |P09{ cannot be treated 



analytically within the proposed scheme and are discussed elsewhere |EHKKLS11| . 

The space-time metric with the choice of parameters (|VIII.8|) considered here is possibly not astrophysically or 
cosmologically relevant due to the negative sign of the cos mological constant. But it could be interesting in the 
framework of the AdS / CFT correspondence |M98l . |GKP98 !| . Our motivation to study the motion of test particles 
in this space-time is more of mathematical character. Furthermore, our results concerning the particle motion in 
Hofava-Lifshitz space-times exhibit the same mathematical structure as for a number of space-times mentioned in the 
introduction. 

Using the substitution q — y/r we find that radial part of the geodesic equation is of the form 



1 dg' ^ 



,9 df 
where k = 8 with 



Pk{q) , (VIII.9) 



Psiq) = ^ (eAwq'' + ea{-K^f'S'' + (E^ - e + KwLl)q^ + Lla{-Kwf/\ - hi) (VIII.IO) 
for massive particles, while fc = 4 for massless particles with 

PM = {{E' + hwLl)q' + Lla{-KwY'\ - L^) . (VIII.ll) 

We then find that 

V-Vo= f /J,, . (VIII.12) 
B. Light rays 

For light we have e = 0. We write the 4th order polynomial (jVIII.llj) as P4:{q) — b^q^ + biq + bo. Introducing a 
new variable x with 

q^- + q4 (VIII.13) 

X 

where q^ is any root of P4{y) we find that (IVIII.12|) reduces to 

I r dx 1 r dx ,..TTT..^ 

(p-ipo = — / + - / , (VIII.14) 

<}iJxo VWx) qiJxo {l + qAx)y/]Mx) 

with 

Pz{x) = (61 + 464(j|)x^ -I- Gb^qlx'^ + ib^q^x + 64 =: a^x^ + a2X^ + aix + gq . (VIII. 15) 
Using the substitutions 

x^^z + P, 7={P, ^ = -T^ (VIII.16) 
V 0.3 3a3 
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this can be brought to the Weierstrafi form 



.2L 

94 



dz 



dz 



with 



92 



yj^z'^ - giz - 53 JzQ (1 + 94(72 + /3))^4z3 -522 + 93 



(VIII.17) 



3,4/ 3aia3 — flj 



.93 = -oo + 



dz' 



fla \ 3a3 / ' 803 27a^ 

In order to invert the ehiptic integrals we introduce v such that v — vi^ = J^^ 
V — v — vq— f °° , '^^ Using p'iv) — ^jAp^{v) — g2p{v) — 33 equation (|VIIL17I) can be rewritten in the form 



(VIII. 18) 
Then z — p{v), where 



ip-ipo 



J_ 

94 



dv' 



1 



dv' 



947 7t,o p(^'') - p{^v) 



which becomes 



^p-^PQ^ 

94 



w — Wo + 



1 1 

947 p'(l'p) 



2C(«p)('^-«o) + ln 



(j{v + i;p) cr(uo + Vp) 



(VIII. 19) 



(VIII.20) 



where is defined by p(fp) = — i±2i^. The solution (jVIII.20[) gives ip = ip{v) and the inversion yields v = v{ip). We 
can then find z{ip) and, thus, r(ip), by substituting v into z = p(''^)- 



C. Motion of massive particles 

For massive particles we have e = 1 and, thus, k = 8. Introducing the new coordinate z through 

1 

9 = - + 98 , 
z 



where qs is any root of ^3(9) we find from (|VIII.12p 



z^dz 



zdz 



dz 



qjzo y/Priz) 98 Jzo y/P7{z) 9| J zo it J zo {Qs' + z)^M^ 



dz 



The curve w'^ — Pjiz) is a hyperelliptic curve of genus g — i. We then introduce v such that v — vq ^ J. 



The solution of this integral is z{v) 



(see Eq. dULlD), where 



V~Vq 

u = ^,+ \ .h{v~vo) 

f2{v - Vq) 



/l(0) = /2(0) = 



(VIII.21) 
(VIII.22) 

2 dz' 
^0 ^Pr{z')' 

(VIII.23) 



and where the functions fi{v ~ vq) and /2(f — fo) can be found from the conditions (7{u) = and (7:^( u) = 0. Als o 
Zo is chosen as a branch point of the polynomial V7{z) which defines the half-integer characteristic 2ti |EHKKL"ii| . 
Through a comparison with the Ui = J dui from pi.ll) we obtain from (|VIIL22p 



1 1 1 I 

ip-Lpa = h{v - vq) + —fi{v - Wo) 3 (w - Wo) + — / 7— T 

98 98 98 % Jzo (qs 



dz 



{qs'+z)^IM7) 



(VIII.24) 



Here the last differential in the equation above is of the third kind and was discussed in section IIII CI for arbitrary 
genus of the underlying polynomial curve (see equation (jllLlTp V Here Z — —q^^, W — \J P-j{Z^. Then the solution 
of (IVIH.22I) is 

<^ - Vo = - — !i{v - Wo) + ^/i(« - wq) - i(w - Wo) 

98 98 98 
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qtW 




(ZM) 



-In 



aiW-iz)) 



In 



(62,0) 



Ve)--3iZ,W) 



a{W+{z)) c7{W+{zo)) 



(VIII.25) 



where W±{z) = J^du± /(^^'J^^ du - K^o 

The solution (|VIIL25p represents a generahzation of the case of genus 1 given in (jVIII.20[) . 



IX. CONCLUSION AND OUTLOOK 



In this paper we developed the inversion of general hyperelliptic integrals of the first, second and third kind. Besides 
that, in pi. 601) we explicitly solved the integration of meromorphic differentials in terms of the ^-function, the 3~ 
vector and the half-periods tj and r)' . Moreover, we provided a method which reduces the number of periods which 
need to be calculated explicitly. We pointed out that computer algebra should be used for the calculation of the 
period matrices in any arbitrary basis which provides a quick and convenient method for the calculation of the matrix 
K and the meromorphic half-periods. For this one needs to know the components of the matrix p,;j which can be 
easily calculated with the help of Lemma IIV.2I 

As a first example we applied this method for solving the geodesic equation in particular cases of Hofava-Lifshitz 
black hole space-times. We considered special cases related to underlying algebraic curves of genus one and three 
and presented the associated analytical solutions for the geodesic equations of massless and massive test particles. 
Other examples where this m ethod will be applied are geodesies in Myers-Perry space-times MP86{ and in black 
ring space-times |ER02l iPSOej . 

When trying to solve the geodesic equation for general Hofava-Lifshitz black hole space-times, that is, for C2 7^ 
in (IVIII.7[) . there is no way to get rid of the square root. In this case one has to square the whole equation, thus 
arriving at a differential equation which is based on a quartic algebraic curve. Similarly, quartic problems also appear , 
for instance, for the geodesic motion in string theory inspired space-times such as Gaufi-Bonnet space-times [BD85j |. 
as well as for t he motion of charged particles in the space-time of the regular black hole given by Ayon-Beato and 
Garcia jABQa iGHKKLMll| . 
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